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ABSTRACT 

We investigate the ultraviolet and optical properties and environment of low redshift 
galaxies detected in the Herschel Astrophysical Terahertz Large Area Survey (H-ATLAS) 
science demonstration data. We use the Sloan Digital Sky Survey seventh release and the 
Galaxy And Mass Assembly database to select galaxies with r Pctro ^ 19.0 mag in the redshift 
range 0.02 ^ z ^ 0.2 and look for their submillimeter counterparts in H-ATLAS. Our results 
show that at low redshift, H-ATLAS detects mainly blue/star-forming galaxies with a minor 
contribution from red systems which are highly obscured by dust. In addition we find that 
the colour of a galaxy rather than the local density of its environment determines whether it 
is detectable by H-ATLAS. The average dust temperature of galaxies that are simultaneously 
detected bv both PACS and SPIRE is 25 ± 4 K, independent of environment. This analysis 
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1 INTRODUCTION 

The wide range of observed physical characteristics of galax- 
ies is indicative of a large variance in galaxy formation his- 
tory. Observations show that the environment of a galax y plays 
an im portant role in shaping its observed properties. iDresslerl 
(1980) showed that galaxy morphology is a strong function of 
galaxy density and numerous studies since then have demon- 
strate d the dependence of galaxy proper t ies on local enviro n- 
ment dLewis etal. 2002; Gomez et alj2003; Balogh et al. 2004315 



Blantonetal. 2005; Park etal. 2007; O'Mill, Padilla & Lambas 



20081; iBall. Lovedav & Brunneij2008l ; lLee et alj201fj) . 

Results from large sky surveys such as the Sloan Digital Sky 
Survey (SDSS) have revealed that the distribution of galaxy colours 
is bimodal, i.e., the so called 'blue cloud ' versus the 'red sequence' 
dStrateva et al]|200ll ; iBaldrv et alj|2004) . The relative numbers of 
blue and red galaxies at a fixed luminosity are observed to vary 
strongly with local density, with blue and red galaxies predomi- 
nantly populating low and high density environments, re spectively 
(e.g. lBalogh et al.l2004bl ; [Ball, Lovedav & Brunnerl2008l) . The ma- 
jor factors that influence the observed colour of a galaxy are its star 
format ion history (SFH), the amo unt of dust attenuation, and metal- 
licity (Johnson et al. 2006, 2007). The vast majority of galaxies in 
the blue cloud are actively forming stars while the red sequence 
consists mainly of early-type passive galaxies (i.e., having little or 
no current star formation) with additional minor contributions from 
heavily obscured star-forming galaxies or edge-on systems. 

Analyses of the dust attenuation in star-forming galaxies sug- 
gest that in comparison to quiescent systems, star- forming ob- 
jects are heavily affected by internal dust extinction dWvder et alj 
20071; Ijohnson et alj|2007l ; [Driver et alj|2007l ; ICortese et alj|2008t 
Toieiro et alj |2009|) . In such galaxies, UV radiation heats dust 



grains which then re-radiate at far-infrared (FIR) wavelengths. 
However, thermally emitting dust detected by, for example, IRAS 
(at 25 — 100/im) can often only constitute a small fraction of a 
galaxy's total du st mass. Results from t h e SCUBA Local Universe 
Galax y Survey dDunne & Ealesl 1200 ll ; Iviahakis. Dunne & Ealesl 
20051), as well as tho se of the sample of the V irgo cluster studied in 



Popes cu et all d2002h and iTuffs et all d2002h . using the ISOPHOT 



instrument on board the Infrared Space Observatory (ISO), show 
that there is a population of galaxies containing much larger pro- 
portions of cold dust that radiates at > 100 pm in submillime- 
ter bands. As such, submillimeter observations provide invaluable 
information for estimating the total dust content of these galaxies 
which in turn helps us to understand relationships with other galaxy 
properties. For instance, although we know there is a clear division 
between the optical properties of galaxies on the red sequence and 
those in the blue cloud, we do not know if this is reflected in their 
submillimeter properties. 

In the present study, we use submillimeter data acquired by 
the Herschel Space Observatory as a part of the H erschel As- 
troph ysical Terahertz Large Area Survey (H-ATLAS; lEales et alj 
l2010l) . The H-ATLAS survey is an open-time programme to sur- 
vey ~600 deg 2 of the extragalactic sky over the wavelength range 
110-500 /im u sing the Photodetecto r Array Camera and Spec- 
trometer (PACS; Pog litsch etal]|2006ri and Spectr al and Photomet- 
ric Imaging REceiver (SPIRE; iGriffin eTai]|2010ri . H-ATLAS also 
covers the largest area in a blind extragalactic survey with Herschel. 
One of the aim if to provide the FIR equivalent for t he SDSS and 
2dFG RS (Two-degree Field Galaxy Redshift Survey; Colle ss et alj 
l200ll) surveys is to characterize the properties of nearby galaxies. 
Also, for all regions covered by H-ATLAS, spectroscopic com- 



pleteness will be guara nteed only up to the redshift of the GAMA 
survey Q(e.g. z < 0.2; iDriver et al.ll2010h . So detailed statistical 
analysis will mainly be possible in the local universe. Therefore it 
is crucial to understand what kind of nearby galaxies are detected 
by H-ATLAS to provide a useful characterization of the FIR prop- 
erties of local galaxies. In addition, we are also going to charac- 
terize the relationship between the nature and environment of sub- 
millimeter sources detected in H-ATLAS. So by concentrating on 
sources at low redshifts, we are able to have a better estimation of 
the galaxy projected density since such measurements become dif- 
ficult at higher redshifts due to a decrease in the observed limiting 
luminosity of galaxies with spectroscopic redshifts. Also, nearby 
sources can be spatially resolved far more easily (e.g. less blend- 
ing), so that we can see where the submillimeter emission is com- 
ing from, unlike for sources at high redshifts. 

The submillimeter data analysed in this study were acquired 
during Herschel's Science Demonstration Phase (SDP) and cover 
an area of ~ 14.5 deg 2 . Three main advantages of Herschel over 
previous infrared observatories are its broad wavelength range, 
high sensitivity and high angular resolution. The wavelength range 
in particular encompasses the peak of the far-IR/submillimeter part 
of the spectral energy distribution for low redshift galaxies thus al- 
lowing accurate measurements of total dust mass. In contrast, other 
submillimeter surveys such as those conducted by the Balloon- 
borne Large-Aperture Submillimeter Telescope (BLAST) with ob- 
servations at 250/im, 350/im, and 500/im, have lower angular res- 
olutio n and do not cover more than ~20 d eg 2 of the extragalactic 
sky dPascale et al.ll2008l ; lDevlin et alj|2009h . 

The layout of this paper is as follows. In §2, we provide a 
summary of our submillimeter sample selection, optical photomet- 
ric and spectroscopic observations from the SDSS/GA MA surveys, 
and U V observations from GAMAy 'GALEX survey dSeibert et alj 
120101) . A summary of the parameters estimated from our data is 
given in §3. Our final results and concluding remarks are given in 
§4. 



2 DATA 

We will use the SDSS DR7 dAbazaiian etafll2009h together with 
the Galaxy And M ass Assembly redshift survey (GAMA survey; 
IDriver et to compile an optically selected sample of low- 

redshift galaxies, and from there, look at the optical properties of 
those sources detected by Herschel. Here, we summarize the data 
in each survey and describe our sample selection. 



2.1 Optical/UV data 

Optical sour ces used in this paper were initially taken from the 
SPSS PR7 dAbazaiian et alj 120091) . Our aim is to have an op- 
tically selected sample of galaxies with the criterion 0.02 ^ 
z < 0.2. To do this, all SPSS sources from the GALAXY 
table were selected and their positions were cross-matched 
within 0.5 arcsec with objects from the GAMA survey in or- 
der to extract their spectroscopic redshifts as well as colours. 
The GAMA survey is based on data from SDSS DR6 and 
UKIDSS-LAS DR4 and consists of three equatorial regions named 
GAMA-09h, GAMA-12h, and GAMA-15h. The field observed 
by Herschel, coincides with the GAMA-09h equatorial region 



A detailed description of the GAMA survey is given in Sec l2. 1 1 
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(Ra = 129.0, 141.0 deg; Dec = -1.0, +3.0 deg) where its spec- 
troscopic survey is com plete up to Petrosian magnitude, r < 19.4 
( Robo tham et alJ l201oh . Moreover, the GAMA data contains r- 
band defined aperture matched photometry in FUV/NUV (far- 
ultraviolet/ ne ar-ultraviolet) fro m GALEX and in ugriz optical bands 
from SDSS dHill et al.l koilh . We only select objects with the 
GAMA spectroscopic redshift quality para meter Q ^ 3, which ap- 
proves their inclusion in scientific analysis jDriver et alJ|201Ch . In 
our optical sample, ~1.4 per cent of sources with 0.02 ^ z ^ 0.2 
have Q < 3. 

We find a total of ~ 3370 galaxies in the redshift range 0.02 ^ 
z ^ 0.2 within wl2.5 deg 2 of our «14.5 deg 2 SDP field. The 
coverage is not complete, because the GAMA data goes down to 
Dec = — 1.0 deg. In addition to opti cal fluxes, NUV magnitudes 
were available for ~ 2680 galaxies JSeibert et alj|201(it iHill et all 
I20H . All magnitudes are corrected for Galactic extinction. 



2.2 Submillimeter data 

The H-ATLAS SDP data were acquired in November 2009 using 
Herschel's parallel/fast scan mode, scanning at 60 arcsec s~ . In 
this mode, mapping by PACS at 100 and 160 fim and SPIRE at 
250, 350, and 500/im were carried out simultaneously. The SDP 
data covers an area of ~ 4.0 deg x 3.6 deg centered on (a, S) ~ 
(09 h 05 m ,+0°30'). 

A complete description of the data red uction and noise anal- 
ysis of PACS and S PIRE data are given in llbar et ail d201fjh and 
IPascale et all feOlfjh . respectively. Observed time-line data from 
both instruments were processed by using the Herschel Interactive 
Processing Environment (HIPE) using custom reduction scripts. 
Thermal drift in bolometer arrays was corrected by applying high- 
pass filtering. The naive map-making method of HIPE was used to 
project the two cross-scan timeline observations. 

The estimated point spread functions (PSFs) of the final maps 
have full width at half maximum (FWHM) values of 8.7, 13.1, 18.1, 
25.2, and 36.6 arcsec at 100, 160, 250, 350 and 500 (im, respec- 
tively. The estimated la noise levels at 100, 160, 250, 350, and 
500Atm are 25-30, 33-48, 6.6 , 7.6, and 9.0 mjy/beam, respectively 
dlbar et alfeOlfilRigbv et alj20ld) . 

Source s were extracted using the 250 (im map as described in 
iRigbv et all fcoiOh . The typical positional error for a ^5 a source 
is 2.5 arcsec or less. For each 250 fim source, corresponding 350 
and 500 flux densities were measured by using the noise- 
weighted/beam-convolved 350 and 500 /an maps at the source po- 
sition determined from the 250 /im map. Fluxes at 100 and 160 
pim were estimated by matching each 250 (im source to the nearest 
PACS sources within a radius of 10 arcsec. 



2.2.7 The submillimeter source catalogue 

The H-ATLAS standard source catalogue of Ri gbv et ail d2010l) is 
supplemented with cross-identification inf ormation from GA MA 
and SDSS DR7 surveys as described in ISmith et alJ (20ld) . A 
likelihood- ratio analys is i Sutherland & Saunders! 1992 ) is per- 
formed by ISmith et alJ feOlCh to match 250^tm sources to SDSS 
DR7 sources brighter than r=22.4 mag. The probability that an op- 
tical source is in fact associated with the submillimeter source has 
been used to defin e the Reliability parameter. According to 
ISmithetal.ld2010h . objects with Reliability^ 0.8 are consid- 
ered to be true matches to submillimeter sources. 

While we use the SDSS galaxy IDs to obtain the H-ATLAS 



PACS and SPIRE subm illimeter fluxes, we apply the same reliabil- 
ity cut as suggested bv lSmithetal] f2010). i.e. Reliability^ 
0.8 which guarantees that we are ~93 per cent (~ 83 per cent) 
complete in the redshift range 0.02 < z < 0.1 (0.1 < z «S 0.2). 
As such, throughout this paper, we refer to a 'detected source' as 
an SDSS galaxy with a 5. Oct submillimeter counterpart that 
has Reliability^ 0.8. Accordingly, we detect 496 galaxies 
at 250/im of which 482 (~ 97.2 per cent) have associated NUV 
fluxes. Thus in 14 galaxies (< 2.8 per cent) with submillimeter de- 
tections we found no UV counterpart, but this will not affect the 
main results of the paper. 



3 ANALYSIS 

3.1 Summary of source detection 

The first step of our analysis is to determine the rate of source de- 
tection in H-ATLAS as a function of galaxy r-band magnitude. The 
r-band magnitude is useful since it can be used as a starting point 
in sample selection for other extragalactic studies of low redshift 
galaxies using the Herschel submillimeter data. To do so, for all 
galaxies, we plot (u — r) vs. Petrosian magnitude r Pctro corrected 
for Galactic extinction^. 

As Fig.[Tjshows, the overall distribution of galaxies (gray con- 
tour lines) in u — r colour is bimodal, showing blue cloud and 
red sequence galaxies. The overlaid data points represent galax- 
ies with submillimeter detections made by SPIRE (left panel) and 
PACS (right panel). Horizontal and vertical histograms in each 
panel show the distribution of data points corresponding to indi- 
vidual SPIRE/PACS bands along the magnitude and colour axes, 
respectively. 

Fig. [U shows that the majority of galaxies with submillime- 
ter detections are located in the blue cloud, with histograms of 
their colour peak around it — r~1.4 without substantial red tails. 
Thus, within the redshift slice of the current sample of galaxies, the 
survey on the whole is mapping the blue sequence. Furthermore, 
a comparison between histograms at 250 /im and those in other 
SPIRE/PACS bands show that the SPIRE 250 fim band is by far 
the most sensitive, since the estimated rate of detections are higher 
at 250^m compared to other submillimeter bands. This is due to 
the higher sensitivity of SPIRE over PACS as well as the higher 
resolution of the 250^im map. With this in mind, in the analysis 
that follows, we only consider sources with 250 fim detections. 

The quantities plotted in Fig.[Tjare also shown in Fig. [2] but as 
a function of absolute r-band magnitude M r and for sources with 
250 fim detections. 'Left' and 'right' panels correspond to redshift 
bins 0.02 ^ z ^ 0.1 and 0.1 ^ z ^ 0.2, respectively. The absolute 
magnitude M r is given by 



Mr 



5 log 



Dl 
lOpc 



(1) 



where Dl is the lu minosity distance, and k r is the k-co r rectio n 
using the method of Chilingarian, Me lchior & Zolotukhinl (2010). 



2 The SDSS Petrosian magnitude is a modified form of the iPetrosian] 
1 1976) system in which galaxy fluxes are measured within a circular aper- 
ture Tip (Petrosian radius) such that the ratio of the local surface brightness 
in an annulus at IZp to the mean surface brightness within IZp, is equal 
to some constant value. Alternatively we can use the SDSS Model magni- 
tude r M . , though it does not change the results since from a sample of 
10000 galaxies within the magnitude/redshift range of our sample we find 
rp„„„ ~ r\, . , +0.07. 
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The adopted cosmology to estimate Dt, is a flat Universe with 
the matter density parameter f2M=0.3 and cosmological constant 
J7a=0.7, with Ho=70 km s _1 Mpc -1 . The 'horizontal green his- 
tograms' in Fig. [2] show that the rate of detections of galaxies in 
the r-band become incomplete for sources with M r > —21.0 mag 
and M r > -21.5 mag in the redshift bins 0.02 < 2 < 0.1 and 
0.1 ^ z ^ 0.2, respectively with almost double the number of 
detections in the higher redshift slice. Thus in the low-redshift Uni- 
verse and within the current depth of H- ATLAS, most of the Her- 
schel detections in the r-band consist of more luminous galaxies 
(and therefore those with larger stellar masses) which have lower 
values of it — r. 

In order to analyze our sample in more detail, we divide 
galaxies into two populations by applying a cut on their u — r 
colour indices. To find an optimum colour cut, a double-Gaussian 
function is fitted to the histogram for colour of all galaxies with 
r potro ~ 19-0 ma g- The combined function has a minimum around 
u — r w 2.2 mag. As such, from here on, unless otherwise stated, 
we refer to sources with u — r ^ 2.2 mag and u — r > 2.2 mag as 
blue and red objects, respectively. 

In Fig. [3] we examine the differential fraction of detections of 
galaxies at 250^m as a function of r-band magnitude. We do this 
using y/n counting errors (n is the number of detected sources in 
each magnitude bin) for all optical sources (black histogram) as 
well as blue (blue dashed histogram) and red (red dash-dotted his- 
togram) galaxies. The results, as presented in Fig. [5] show that (i) 
In general > 50 per cent of SDSS galaxies with r Pctro ^ 17.0 mag 
have 250/tm counterparts in H-ATLAS. (ii) The submillimeter de- 
tection rate falls off rapidly as one moves toward fainter galaxies. 
This is mainly due to a decrease in the observed signal-to-noise ra- 
tio toward fainter galaxies, (iii) In each bin of r-band magnitude 
and for r Pctro < 18.0 mag, the detection fraction is > 4 times 
higher for 'blue' galaxies when compared to the 'red' population, 
indicating that the vast majority of detected sources are blue ob- 
jects. Although such blue sources are in general star-forming galax- 
ies, it does not mean that 250/im red detected objects are passive 
galaxies. We will discuss the properties of such sources in Sec- 
tion|3~3l 

Finally in Table. [TJ we present the total percentage of sources 
detected by H-ATLAS in samples of SDSS galaxies with differ- 
ing limiting magnitudes. Samples of SDSS galaxies in the second 
column include galaxies with all colour indices, while those in the 
third and fourth columns contain blue and red objects, respectively. 



3.2 Environment of detected sources 

Our results in i)3.1| show that submillimeter detected sources are 
preferentially blue sequence galaxies. Thus, since blue sequence 
galaxies are much more frequent in low density environments, we 
expect to find a lower detection rate in group/cluster environments 
than in the field. 

To explore the environmental density of submillimeter de- 
tected galaxies, we consider the projected surface density, based on 
counting the number of nearest neighbours, i.e. the density within 
the distance to the Nth nearest neighbour. Although this is a 2D 
estimate, the redshift information of each galaxy is used to remove 
the background and foreground sources. The environmental density 
Ejv around each galaxy detected at 250/im is determined from the 
following relation: 

E^Mpc" 2 ) = -*L (2) 
nd N 



where djv is the projected comoving distance to the N 
nearest neighbour t hat is within the allowed redshift range 
±Azc=1000 km s" 1 l lBalogh et alj2004aklBaldrv et alj|2006h . We 
set N = 5 in our analysis. The neighbouring galaxies are those 
with spectroscopic redshifts and M r ^ M t Umit — Q(z — Zq), 
where Q=l.6 is the evolution of th e galaxy r-band lum inosity 
relative to zo = 0.1, determined by Blanton et al. (2003). Hori- 
zontal histograms (grey thick lines) in the 'left' and 'right' pan- 
els of Fig. [2] show that galaxy counts become incomplete for the 
M r > -19.0 mag range if 0.02 sC 2 < 0.1 and M r >-20.5 mag 
if 0.1 ^ 2 ^ 0.2. Thus the parameters {M limit ,2o) were set to 
(-20.5,0.1). 

The rate of detection of the H-ATLAS sources as a function 
of the density parameter E5 is plotted in panel 'a' of Fig. [4] Error 
bars are Poisson errors on the mean. Detection rates have been es- 
timated for galaxies with r Potro 17.0 mag, since galaxies with 
larger r-band magnitudes are not systematically detected due to 
their low signal-to-noise ratios. In other words, the majority of 
fainter galaxies are not detected in H-ATLAS, irrespective of their 
environments. 

It is evident that the fraction of detected sources in low-density 
environments is larger by a factor of two in comparison to the rate 
of detection in high-density environments. This is expected, since 
the majority of detected objects in H-ATLAS are blue galaxies 
that preferentially populate the low-density regions. Results remain 
more or less the same if we use E4 or Eg instead of E5, showing 
that the dependence of the detection rate on environmental density 
is not sensitive to the choice of N in Eq.|2] To obtain a better under- 
standing of the type of detected objects that populate either low or 
high density environments, we plot the rate of detection of galax- 
ies as a function of r Potro for sources with different colours and 
densities in Fig.[4](see panels 'b' and 'c'). 

The estimated projected density ranges from —2.0 to 2.0 on 
a logarithmic scale. We therefore define the threshold between low 
and high density environments at log E5 ~ 0. As such, in pan- 
els 'b' and 'c' of Fig. [4] black histograms correspond to sources 
in low density environments with log(Es/Mpc 2 ) < while red 
dashed histograms refer to the fraction of detected galaxies in high 
density regions, i.e. log(Es/Mpc 2 ) ^ 0. In addition, sources in 
panel 'b' have index u — r ^ 2.2 mag (e.g. blue objects) and the 
colour of those in panel 'c' is u — r > 2.2 mag (e.g. red objects). 
Although the number of blue galaxies is smaller in high-density re- 
gions, panel 'b' shows that the rate of detection remains the same 
in either environment. The situation is similar for the red galaxies, 
although the overall rate of detection of red objects is much smaller 
than for blue ones. 

3.3 Nature of RED 250/im sources 

We have already seen that the majority (~ 80 per cent) of galax- 
ies detected in H-ATLAS at low redshift are blue galaxies while 
the rest ~ 20 per cent are red. However, such a fraction on red de- 
tected sources is an upper limit since many of these galaxies might 
have traces of recent star formation activities which is not reflected 
in their optical u — r colour. In fact, UV-optical colours are more 
robust in classifying galaxies into blue/red categories in compari- 
son to optical colours s ince they are more sensitive to recent star- 
forma tion activity (e.g. ISchawinski et al.ll2007h ICortese & Hughes 
2009). The reason that we did not use UV-optical colour in Sec. 3. 1 
is that for ~ 20 per cent of SDSS galaxies in the main sample, UV 
fluxes were not available from GAMA (see Sec. 2.1). These galax- 
ies (e.g. those without UV counterpart) are basically faint while 
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Figure 1. Colour-magnitude diagram of optically selected galaxies (grey contour lines) and their submillimeter counterparts with Reliability^. 8 de- 
tected at > 5cr by SPIRE (left panels) and PACS (right panel). Horizontal and vertical histograms show the distribution of galaxy optical colour (u — r ) and 
r-band magnitude (r p t ), respectively. Colours represent: 'grey' all galaxies; 'black' sources with SPIRE 250/tm (left panel) and PACS lOO^m (right panel) 
detections; 'pink' sources with SPIRE 350/im (left panel) and PACS 160^im (right panel) detections; 'green' sources with SPIRE 500/^m detections. 
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Figure 2. Same as in Fig.[T]but as a function of absolute r-band magnitude M r and for 250/im detections. 'Left' and 'right' panels correspond to redshift bins 
0.02 ^ z ^ 0.1 and 0.1 ^ z ^ 0.2, while horizontal histograms show the distribution of galaxy u — r colour and absolute r-band magnitude, respectively. 
Colour of data-points/histograms represent: 'green' all galaxies; 'black' sources with SPIRE 250/^m detections. 
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Figure 3. Histograms of the fraction of detected galaxies as a function of r-band magnitude. The black histogram includes all galaxies, while the two other 
histograms represent the detection rate of blue (blue dashed histogram; u — r ^2.2) and red (red dash-dotted histogram; u — r >2.2) sources, respectively. 



Table 1. The total fraction of sources detected by H-ATLAS among SDSS galaxies to differing limiting magnitudes, as shown in the first column. Samples of 
SDSS galaxies in the second column include galaxies with all colour indices, while those of third and forth columns contain blue and red objects, respectively. 



sample limiting magnitude Total fraction detected (percent) 

all colours blue (u — r ^ 2.2 mag) red(ii — r > 2.2 mag) 

r 



Petr 
Petr 
Petr 
Petr 



^ 16.0 mag 53 ± 8 97 ± 18 27 ± 8 

sS 17.0 mag 42 ± 4 75 ± 7 17 ± 3 

sS 18.0 mag 29 ± 2 47 ± 3 11 ± 2 

O9.0 mag 17 ±1 23 ±2 9±1 
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Figure 4. (a) The fraction of SDSS galaxies detected in H-ATLAS as a function of galaxy local density £5 estimated using Eq. [2] (b) Histograms of the 
fraction of detected blue (u — r ^2.2 mag) galaxies as a function of r-band magnitude. 'The red dashed histogram' includes galaxies in high density regions 
(log(S5 /Mpc 2 ) ^ 0) and 'the black histogram' includes those in low density regions (log(Ss /Mpc 2 ) < 0). (c) Same as in 'b' but for red (it — r >2.2 mag) 
galaxies. 
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Figure 5. (a) NUV— r colour distributions for all the SDSS galaxies (black histogram) with available UV fluxes as well as those detected in 250/im (green- 
shaded histogram). (b)The NUV— r vs. vL25o/^^nuv f° r galaxies with 250/^m detections (green triangles). Overlaid (black diamonds) are galaxies detected 
in all SPIRE bands. 



their optical colour distribution peaks around u — r ~ 2.5. So in or- 
der to avoid any biases in analyzing the fraction of source detection 
in previous sections, we were restricted to use the SDSS optical 
colours. Fortunately, NUV fluxes are available for > 98 per cent 
of detected sources in 250/tm which enable us to investigate the 
spectral properties of galaxies with 250/mi detections from their 
UV-optical colours. 

Thus in Fig. [5^, the NUV— r colour distributions for all the 
SDSS galaxies with available UV fluxes (black histogram) as well 
as those detected in 250/im (green-shaded histogram) are shown. 
The black histogram in Fig.|5^ shows a bimodality with a minimum 
around NUV— r w 4.5 mag found from fitting a double-Gaussian 
function to it. In addition, ~ 95 per cent of detected sources have 
NUV— r ^ 4.5 mag (i.e., blue in UV-optical colour) and just a 
small population of detected sources have NUV— r ^ 4.5 mag (i.e., 
red in UV-optical colour). 

It is therefore interesting to know whether such red objects 
(i.e., NUV— r ^ 4.5 mag) are obscured star-forming systems or 
whether they are a population of quiescent/passive galaxies. As 
such, for all detected galaxies, a UV-optical vs. infrared colour- 
colour diagram is plotted in Fig. |5p, where green triangles repre- 
sent sources detected at 250/im and highlighted data points (black 
diamonds) represent those H-ATLAS sources detected in all three 
SPIRE bands. Not surprisingly, the plot shows a correlation be- 
tween NUV— r and vL25o/vLnuv, with red galaxies having a 
higher ^I/25q /^Lnuv ratio t han b lue systems. A similar relation 
was found by I Johnson et al.l d2007l) between the total far-infrared- 
to-near-ultraviolet luminosity ratio (I/tir/^nuv) and NUV— r for 
a sample of local galaxies selected from SDSS. 

By combining this re lation with the l i brary of infrared spec- 
tral energy distributions o f ICharv & Elba3 fcOOll) . we can investi- 
gate the real nature of the H-ATLAS sources with NUV— r >4.5. 



According to ICharv & Elbaj d200ll) . the ^Z/25o/£tir luminosity 
ratio varies in the range ~4-65. Even assuming the lowest pos- 
sible value (i.e., ~4), this implies a Ltir/^nuv ratio ~6-250 
for our red galaxies and we can use this value to estimate their 
NUV dust attenuation A(NUV). If we adopt the typical rela- 
tion between I/tir/^ nuv and A(NUV) for star-forming objects 
(e.g. lBuat et a i]|2005h . we find A(NUV) ~l-4 mag. Interestingly, 
even usin g the relations betwe en LtirALnuv and A(NUV) pre- 
sented bv lCortese et alj ( 120080 which take into account the heating 
of dust coming from evolved stellar populations, we still obtain 
A(NUV) ^0.2-l.fl 

Thus, this simple exercise shows that not surprisingly, the vast 
majority (^90 per cent) of the red objects detected by H-ATLAS 
in our sample have 'corrected' NUV— r colour lower than 4.5 mag, 
i.e., they are red not because they are old/passive but mainly be- 
cause they are obscured by dust. Of course, the presence (even if 
low in fraction) of a population of 'truly passive' systems in the 
H-ATLAS survey is a very intriguing possibility and might have 
important implications for our understanding of dust properties in 
galaxies. A detailed search for FIR emitting 'passive' systems and 
a discussion of their p roperties will be presented in a future work 
jRowlands etalj201lh . 



3.4 Variation of temperature with density 

Our findings in Section [T2l show that, irrespective of their environ- 
ments, the rate of detection of submillimeter sources remains the 



3 These values have been ob tained by using the most extreme case dis- 
cussed bv lCortese et alj (2008), i.e. r < 2.8, and they are thus lower limits 
to the real values for our sample. 
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same in either high or low density regions. Therefore it might be in- 
teresting to test whether the submillimeter properties (e.g., colours 
and dust temperature) of detected sources in H-ATLAS correlate 
with the environment. 

Out of 71 galaxies with ^ 3u detections in all SPIRE bands, 
38 have associated 3<r) PACS detections. To estimate the dust 
temperature of these 38 galaxies, we fit the submillimeter photom- 
etry with a single component modified black-body SED over the 
wavelength range 100-500/im. The equation for a single compo- 
nent black body is given by 

F v = ^MB U {T). (3) 

Here M is the dust mass, T is the dust temperature, B v (T) is the 
Planck function, Dl is the luminosity distance to the galaxy, and 
k v is the dust emissivity. The dust emissivity is assumed to be a 
power law in this spectral range, where k v oc . As we are not 
constraining dust masses, an absolute value of n is not needed. We 
assume /3 = 1.5, though we note that the value is notoriously un- 
certain. 

The uncertainties in the fl ux densities adopt ed from the stan- 
dard H-ATLAS catalogue (e.g. lRigbv e t al. 2010) include the con- 
tribution from instrumental and confusion noise only. Thus we add 
the calibration errors of 10-20 per cent in the 100-160/xm PACS 
bands and 15 per cent in all SPIRE bands in quadrature to those 
given in the H-ATLAS catalogue. The best-fit solution is found by 
minimizing the chi-squared (x 2 ) function, yielding an average un- 
certainty in the fitting of ±1.8 K for temperature. The average re- 
duced x 2 value for all the fits is 0.7, with a standard deviation of 
0.6. We find that the data in the 100-500/im range are accurately 
fitted by a single black-body component. 

Fig. Hjji shows that the submillimeter colour Sinn /S500 cor- 
relates with the dust temperature. Such a correlation also exists 
between the dust temperature and other submillimeter colour in- 
dices, e.g., S160/S350, S250/S500, etc. However, not surprisingly, 
the relation between the S100/S500 A ux rat i° an d temperature is 
the strongest one. The estimated Spearman correlation coefficient 
in panel 'a' is Rs = 0.96 and the best-fit linear regression to the 
data points (black dashed line) suggests the following relation for 
the dust temperature: 

T{K) = 0.51(±0.02) x Sioo/S 5 oo + 19.5(±0.3). (4) 

A plot of the estimated dust temperature as a function of the 
density parameter E5 for all galaxies in our sample is presented in 
panel 'b' of Fig. [6] Galaxies are divided into two subsamples based 
on their stellar mass content with circles and squares corresponding 
to sources with low (9.5 < log(M,/h~ 2 M ) < 10.5) and high 
(10.5 < log(M„//i- 2 M ) < 11.5) stellar masses, respectively. 
Galaxy stellar masses, M*, have b een computed in u nits of solar 
mass from the relationship given by Yang et ah (2007) as follow: 

log(M„//i" 2 M ) = -0.406 + 1.097(c/ - r) 

-0.4(M r - 5 log ft -4.64), ' ^ 

using the (g — r) colour and where M r is the r-band absolute mag- 
nitude corrected for extinction and redshift. The green dashed line 
shows the mean temperature of the whole sample (T» 25 ±4 k]3, 

4 The estimated mean temperature is T mean ss 22 ± 4 K in the case of 
= 2.0. In fact our data show a strong correlation (Rg = 0.99) between 
the estimated temperatures based on either values of /3 such that Tp—2.0 = 
0.78(±0.01) x Tp =l * + 2.1(±0.2). 



while green dotted lines indicate the standard deviation. This result 
fairly agrees with the rest-frame dust tempe ratures measured for a 
sample of the H-ATL AS galaxies studied in lDve at al.l d20ld) and 
lAmblard et alj ( 120101) as well as those of the sample of BLAST 
sources determined bv lDveatal.N2009h . We note that the densities 
in panel 'b' of Fig.|6]span values typical of field galaxies and poor 
groups, whereas we are not able to trace very high density environ- 
ments, such as compact groups or clusters of galaxies. 

From inspection of Fig. [6p, it appears that the dust tempera- 
ture does not show any systematic variation with the local density 
of either low mass or high mass galaxies. Further checks to see 
whether there is any correlation between either dust temperature or 
E5 with other parameters such as galaxy redshift, physical size, and 
luminosity reveal none. 

Interestingly, Fig.|6]shows the presence of six outliers (marked 
as 'A' to 'F') characterized by dust temperatures (T ~ 32 ± 2 K) 
significantly higher than the values typically observed in the rest 
of our sample. SDSS images (see Fig. |7J reveal the presence of 
diffuse tidal stellar features associated with at least three (A, C and 
F) of the six outliers in Fig. [6] These are the only galaxies in our 
sample of 38 objects with obvious signatures of tidal disturbance, 
suggesting that the dust temperature has recently been enhanced 
via gravitational interaction. 



4 CONCLUSIONS 

This study presents a variety of statistics based on Herschel ob- 
servations acquired during the H-ATLAS Science Demonstration 
Phase. We started with a sample of SDSS-GAMA galaxies at low 
redshift (0.02 < 2 < 0.2) and then characterized the properties of 
the matched submillimeter sources detected in the H-ATLAS sur- 
vey from their colours. Our results can be summarized as follows: 

• At least « 50 per cent of all SDSS galaxies with r Pctro < 
17.0 mag are detectable in H-ATLAS. In addition, the vast majority 
of detected sources in H-ATLAS are blue/star-forming objects. The 
majority of undetected faint/blue sources with r Pctro > 17.0 mag 
remain undetected due to their lower signal-to-noise. By using the 
NUV— r index as a proxy for galaxy colour, we detect a small frac- 
tion (~ 5 per cent) of red objects (NUV— r > 4.5 mag) which are 
mainly galaxies with a higher level of dust attenuation. As such the 
red sequence population which is dominated by less-dusty/passive 
objects, do not seem to have a significant contribution to the ob- 
served submillimeter emission in red objects. 

• The rate of detection of galaxies decreases from ~ 70 per cent 
to 30 per cent as one moves from low density to high density re- 
gions. This is due to the fact that low density regions are more pop- 
ulated by blue/dusty objects. However the detection rate of blue/red 
galaxies remains constant in either low or high density environ- 
ments. In other words the colour of an object rather than its local 
density determines whether it is detectable in H-ATLAS. 

• The estimated dust temperature in galaxies with detections in 
all SPIRE bands and > 3a detections in PACS bands is 25 ± 4 K, 
regardless of the environment. However, we show that gravitation- 
ally perturbed systems have temperatures significantly higher than 
the rest of our sample. 

The next crucial step is to extend our analysis to a wider range 
of galaxy densities/stellar-masses to see how far the environment 
of galaxies could affect their observed submillimeter properties. 
While Herschel observations of the Virgo cluster are starting to 
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Figure 6. (a) Dust temperature vs. submillimeter colour 5"ioo / 5*500 ■ The dashed line shows the linear regression fit as described by Eq.|4] (b) Dust temperature 
as a function of the local density of galaxies, E5, for sources with ^ 3cr detections in SPIRE/PACS bands. The green dashed line intersects the y-axis at the 
mean temperature of the whole sample (T ss 25 K) while green dotted lines represent the standard deviation of the mean temperature. Circles and squares 
represent galaxies with low (9.5 ^ log ( M, / ft, ~ 2 Mq ) ^ 10.5) and high (10.5 ^ log(Af*/ft~ 2 Af ) sg 11.5) stellar mass, respectively. The SDSS images 
of objects marked as 'A' to 'F' (with 'A' having the maximum temperature) are shown in Fig. [7] 



J09040 1.02*0 12729.1 .1085748.06+004641.3 J090034. 18*011857.2 J090746.9 1*003430.4 J090734.74+012504.1 J091205.B4+002656.0 





w 

■ 

* . 


1 1 

• 

■ 


# 




• 

• 


A. 


B 


c 


D 


t 

E 


F 



32 kpc 



31 kpc 



24 kpc 



20 kpc 



28 kpc 



21 kpc 



Figure 7. SDSS postage stamp images in the r-band (70 X 70 arcsec 2 ) of galaxies corresponding to the marked data points in Fig. [6] Images are sorted 
(left-to-right) according to the estimated dust temperature. The SDSS optical major isophotal diameter in the r-band (e.g. isoA-r) has been used to measure 
the optical linear diameter of each galaxy D25 in units of kpc as indicated on each panel. 



reveal how t he environment can affect the dust properties of clus- 
ter galaxies dCortese et ai]|2010l) . only a wide-area survey like H- 
ATLAS will make it possible to unveil the effects of nurture across 
the whole range of densities (from voids to the cores of clusters). 
Thus, once completed, the H-ATLAS survey together with data 
from the GAMA survey will offer a unique data set, spanning from 
the UV to the submm, with which to investigate the star formation 
history of galaxies over a large range of galaxy densities/stellar- 
masses. At the same time, properties such as dust temperature and 
dust-mass/stellar-mass of galaxies can be explored as a function 
of environment using PACS/SPIRE submillimeter fluxes. This will 
enable us to further investigate the correlation between star forma- 
tion history and dust properties of galaxies as a function of stellar 
mass and local density. 
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